Phases of a polar spin-1 Bose gas in a magnetic field 
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Bose-Einstein condensed spin-1 gases have rich mag- 
netic properties [l|, 0, H, EL Hi- They can be classified 
according to the sign of the spin-dependent part of the 
interaction. In this paper this sign is assumed to be pos- 
itive, i.e. the interaction is of antiferromagnetic type. In 
the absence of a magnetic field the system at zero temper- 
ature is in the polar state [1, 0] and no spontaneous mag- 
netization is created. An external magnetic field induces 
a magnetization and as shown by Ohmi and Machida [7] 
in the Bogoliubov approximation a complete spin order 
appears when the magnetic field is strong enough. 

The purpose of the present paper is to extend the in- 
vestigations for a homogeneous system to finite tempera- 
tures below the Bose-Einstein condensation (BEC). We 
shall denote by P2 the phase originated from the polar 
state when a small magnetic field is switched on. This 
phase goes over to another one denoted by PI when the 
magnetic field B and/or the temperature is raised. In the 
P2 phase two continuous symmetries are broken while in 
phase PI only one. The ground state spinor has two 
and one components, respectively. Correspondingly, two 
gapless (Goldstone) modes are expected to exist in P2. 
One of them develops a gap when entering the phase PI. 
Extended Hugenholtz-Pines theorems will be presented 
whose fulfilment is a necessary condition for such a be- 
haviour. Another interesting aspect of this phase transi- 
tion is the crossover from a classical to a quantum phase 
transition when the temperature goes to zero, which will 
also be investigated. Furthermore as in [l|, Q a nonzero 
magnetic chemical potential fi m is included in the treat- 
ment. The discussion above remains valid, if the mag- 
netic field B were replaced by B = B + fi m , which can be 
regarded as an effective field determining the magnitude 
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of the average magnetization. 

For demonstration we choose a model which bears 
many features of the full microscopic description. As 
a matter of fact, it has been found in the case of a scalar 
Bose gas that the dielectric formalism can serve as a guide 
to find approximations which meet conservation laws and 
ensure that excitation branches associated with Gold- 
stone modes are gapless @, Q. Such an approach will 
be adopted here in a necessarily extended form since col- 
lective motions now include different types of spin waves 
besides particle density oscillations. Further differences 
arise due to the presence of the effective magnetic field 
B. Different regions are distinguished in the model on 
the B — T plane. Properties are discussed particularly in 
those regions, which might be experimentally accessible 
(note that the theory worked out for a homogeneous sys- 
tem can be applied for a large system as a local density 
approximation) . 

We consider a translationally invariant system of spin- 
1 particles in a box with volume V in a homogeneous, 
external magnetic field pointing to the z-direction. The 
Hamiltonian takes the following form: 

ft = ^ [( 6k ~ t l ) 6rs ~ ^ L ( F z)ra <4(k)a s (k) 

k 

+ ^7 E alik^MV^'aMaAk,), (1) 

k l+ k 2= k 3 + k 4 

where a£(k) and a r (k) are creation and destruction oper- 
ators, respectively, of one-particle plane wave states with 
momentum k and spin projection r. The spin index r 
refers to the eigenvalue of the ^-component of the spin 
operator and can take the values +,0,—. Summation 
over repeated indices is understood throughout the pa- 
per. In this usual basis the spin operators are given by 
F z = diag(l,0,-l) and F x = {F + + FJ)/2 and F y = 
{F+ - F-)/2i, and finally: F± = V2{6, %± 5 S . + 8 r . 5 s , T ). 



In Eq. (p} ek = h 2 k 2 /2M refers to the kinetic energy of 
an atom (M is the mass of an atom) , /i denotes the chem- 
ical potential. The quantities B and cDl are introduced 

as 

g/J, B B = fub-L = fthJL + gilB^rn- (2) 

Here fio>L is the Zeeman energy shift in a magnetic field: 
fiojL = 9^bB, where g is the gyromagnetic ratio, /ib is 
the Bohr magneton; B is the modulus of the homoge- 
neous magnetic field and /x m plays the role of a Lagrange 
multiplier for the magnetization. In Eq. (p} V^ s s is 
the Fourier transform of the two particle interaction po- 
tential, which for the low temperature, dilute gas can be 
modeled by the momentum independent amplitude given 
for spin-1 bosons by 0,(1,01 : 

V?'/ = c n 5 rs S r , s , + c s (F) rs (F) rV , (3) 

In this paper we consider systems with c s > 0. An ex- 
ample of such a system is the gas of 23 Na atoms [Icj . 

It is important that preparing the system with a suit- 
able magnetization the effective field B can be made 
much smaller than the external magnetic field B. This 
procedure made possible to observe experimentally the 
phase P2 (see Fig. 3c in Ref. Q and Figs. 21c, 22c in 
Ref. Q). Note that the quadratic Zeemann shift makes 
the phases PI and P2 unstable for increasing magnetic 
field. It is assumed throughout this paper that the system 
is away from this stability border and the condensation 
occurs only in the spin directions + and — . 

The Hamiltonian (JTJ) is invariant under the gauge 
transformations a± — > a±e lv± and ao — > aoe t 2 t - v + +ip -\ 
This is equivalent to a r — > a r e l ^ +rv \ with <fi = (ip + + 
ip-) I 2 and ip — (ip+—(p-)/2. The invariance under trans- 
formations with if = yields a conservation law for the 
particle number, while with <j> = results in the conser- 
vation of the z-component of total magnetization. 

We investigate such a parameter region (fi,B,T), 
where the system has a Bose-Einstein condensate 
(a r (0)) = \fNcCri with iV c the number of atoms in the 
condensate and ( r = (C-LxCotC-) being the normalized 
spinor of the condensate One can take Co = 0, which 
can be shown to be consistent with the general theory, 
i.e. to all orders in the perturbation expansion. Note 
that when the magnetic field and the magnetization are 
zero the magnitudes of the two components of the spinor 
(C+ and £-) are equal. This spinor is equivalent to what 
is usually taken for the polar state, i.e. C = (0, 1,0) [6j. It 
is convenient to define a new set of creation and annihila- 
tion operators: b r (k) = a r (k) — %/7VcCr^k,o and consider 
the Hamiltonian |T]) as expressed with this new set. The 
canonical transformation, together with the requirement 

(6 r (k)) = (6t(k)> = (4) 

let us to define an ensemble with density matrix p = 
cxp(— f3TL) I Z that exhibits the symmetry breaking asso- 
ciated with Bose-Einstein condensation. 
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FIG. 1: The phase diagram of the antiferromagnetic spin-1 
Bose gas obtained in the mean field approximation 

The general features of the phase diagram can be read 
in Fig. [T] (calculated in an approximation to be intro- 
duced in the second part of the paper). The red line 
designates the border for BEC of the spin-1 gas, T c de- 
noting the critical temperature at B — 0. In phase P2 the 
condensate is characterized by a two component spinor 
C = (C+, 0, C-), IC+| 2 + IC-| 2 =' 1, while in PI C = (1,0, 0). 
For the P2-P1 phase transition the order parameter is 
^//VcC- ■ In the P2 phase both symmetries expressing the 
conservation of the particle number and that of the z- 
component of the magnetization are broken by fixing tp+ 
and if- . In the PI phase like in the ferromagnetic case [6| 
the only remaining phase </? + is related to the breaking of 
the combined symmetry of particle number conservation 
and magnetization. 

The key quantity in our presentation is the finite- 
temperature Green's function of the system: 

^S(k > r)=-(r T [6?(k I r)6f (k,0)]). (5) 

t is the imaginary time and T T refers to the r ordering 
operator pi The Greek indices take the values ±1, 
with 67 (k) = 6 r (k) if 7 = 1 and bj (k) = 6+(-k) if 7 = 
— 1. Here Q{\ (k, r) and Q™x _i(k, r) stand for the normal 
Green's functions, while ^ s _ 1 (k, r) and G™i r) for 
the anomalous ones, which arise due to Bose-Einstein 
condensation. 

The propagator Q® ( ® describes dynamics corresponding 
to rotations around axes lying in the x-y plane. Such 
spin-wave modes have a gap due to the nonzero effec- 
tive magnetic field B. One can convince oneself, that the 
specific choice of the spinor, i.e. Co = 0, and of the di- 
rection of the magnetic field leads to Q™ s = for r = 
(s = 0) and s = ± (r = ±). As a consequence (?™ de- 
couples from the other Green's functions and will not be 
discussed here further since we concentrate on the order 
parameter related to the P1-P2 transition and to its cor- 
relations described by Q~^~ .The structure and couplings 

of the Green's function Q~^~ are different in the phases 
PI and P2. The main results are as follows. 

In P2 g~f is coupled to £+/ by g+f . The 4x4 matrix 
with elements Q™ s , r, s = ± denoted by Q_ takes the form 
in the Matsubara representation 

£ _1 (k,iw„) = £^)(k,iw n ) - E(k,iw„), (6) 
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FIG. 2: (a) The symbolic representation of an irreducible 
diagram having i condensate lines with r = + and j with 
r = — . (b) A fourth order example is shown with i = 4 and 
j = 2. The dashed line represents the interaction, the circle 
stands for the square root of the corresponding condensate 
density and the line with an arrow denotes the free Green's 
function. 



with the self-energy given by the following hypermatrix 
notation: 



£ = 



E++ £ H 

£-+ xr 



1,1 

s 



^1,-1 



1.1 



-1,-1 



(7) 



The inverse of the free propagator is £L Q » 



= diag(S + 
h~ l {e^-pb). 



WL,S*+(Ii L ,S-a) L ,^,*-a) L ), with 5 = iu n - 
The asterisk denotes complex conjugation. 

Concerning the related spectrum it is important that 
the two components of the condensate are independent 
in the sense that they have independent phases: £' = 
(e^+C+,0,e^-C-) is equivalent to ( = (C+,0,C-)- Con- 
sequently two gapless Goldstone modes are expected. In 
this respect it is important that generalized Hugenholtz- 
Pines (GHP) theorems are valid. They are as follows 



£++(o, o) - £+^(0, o) = rr V + sh,, 
sr; il -(o,o)-sr;- 1 (o,o) = ^- 1 /i-w L , 

£+7(0,0)- £+^(0,0) = 0, 



(8b) 
(8c) 



One can show that these relations provide necessary con- 
ditions for the existence of gapless modes, since the de- 
nominator of the Green's function matrix as defined in 
Eq. ([BJ is zero at k = 0, iuj n = if Eqs. JHJ) are fulfilled. 

We sketch the derivation, because it sheds light on 
the structure of the perturbation expansion. Namely, 
all n th order diagrams contributing to the self-energies 
£if (0, 0) and S^ s _ 1 (0, 0) can be obtained easily from the 
one-particle irreducible graphs without external lines but 
with i number of condensate lines with spin projection + 
and j number of condensate lines with spin projection — . 
Particle number and spin conservation results in i and j 
being even. For an illustration of such a graph see Fig. 
[2j where the condensate lines are represented by circles. 
Let us denote its contribution by $a . Here n refers 
to the order in perturbation theory and the index a is 
for distinguishing between graphs with the same n,i,j 
but different topologies and consequently different con- 
tributions. To obtain the n th order contribution of the 



normal self-energy £+-, (0, 0) one has to replace two ap- 
propriately chosen condensate lines by an incoming and 
an outgoing particle line. Note that those circles where 
incoming lines can be substituted are gathered at the left 
side of Fig. [2] (a) while those standing for possible places 
of outgoing lines at the right side. Finally one arrives at 



E++(0,0) = 



1 



EE 



(n)i,j 



(9) 



where n c<r = N c \£ r \ 2 /V is the condensate density in spin 
projection r G {+,—}. For obtaining the anomalous 
£^^(0,0) self-energy one has to replace two appropri- 
ately chosen condensate lines by two outgoing lines. One 
gets that 



i(O,0) 



1 



EE 



- 1 $ 



(10) 



It is possible to calculate this way the contribution of 
the irreducible tadpole diagrams as well, i.e. those with 
one incoming particle line: 



1 



EE 



Using Eqs. ©, {TOD and (HTJ one arrives at 
S++(0,0)-E+± 1 (0,0) = 



(11) 



(12) 



The requirement formulated in Eq. (j4]) for r = + 
satisfied when Eg = ft -1 A 1 + ^l, that yields with |T 
one of the GHP theorems (f5aj) . One can prove the others 
similarly. 

Crossing the border line one enters the PI phase where 
the condensate spinor has only one nonzero component 
C = (1,0,0). This spinor has only one phase, which can 
be chosen freely; a continuous symmetry is restored and 
therefore only one Goldstone mode exists. Furthermore 
in phase PI E+^T = E~+ = 0. Consequently C/++ and 
Gay are no more coupled. The remaining Goldstone 
mode shows up in the spectrum of and the HP the- 
orem £^|(0,0) - £^(0,0) = + w L ensures its 
gapless nature. The order parameter dynamics as de- 
scribed by Gay further simplifies into 

[GifiKiuJn)}^ 1 = iw„-fi -1 (e k -Ai + ^L)-£rf > (I 3 ) 

since S^Zi = fl2| . The excitation spectrum deter- 
mined by Eq. (|13p has a gap, which disappears when 
reaching the phase boundary between PI and P2: 



fiE 



li 



(o,o) L .„. 



(14) 



At T = this condition designates a quantum phase 
transition and it is expected that there is a crossover 
from classical to quantum phase transition when T — > 
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13]. One can define the shift exponent <fi s by the equa- 
tion B C {T) = B c (0) - wT 1 ^- It is useful then to 
introduce the variable x = T/[B — £? c (0)]^ s . Denoting by 
x c its value on the critical line, one can write the leading 
singularities as follows 



Qif (o,o) = 



C(T) 



x - x c 



T > 0, 



(15) 



B-B c {0) 



70 ' 



T = 0. 



One can show that (j) s agrees with the usually defined 
crossover exponent 4> when crossover scaling is valid [in 
this case C(T) ~ T^ ^ in Eq. fity]. One can specify 
a low temperature region where the energy gap is higher 
than the temperature. Here only exponentially small cor- 
rections due to nonzero temperatures are expected. 

In the following we will discuss static and dynamic 
properties within a kind of a mean field approach. We 
start by writing self-consistent equations for the densities 
of the conserved quantities 

n = n C:+ + n c . + n' + + n' + n'_, (16a) 
m = n Ci + — n c . + n, — n'_, (16b) 

with n being the total particle density and m the mag- 
netization density. The density of non-condensed atoms 
in spin projection s S {+, 0, — } is given by the Bose dis- 



tribution n' s 
and 



v 'Ek^. with n 'k, s 



£k,s = ek + c n n - (1 + s(c s m - Swl) 



1 



(17) 



There exists also a relation following from the require- 
ment (6+(0)) = 0, which provides the equation 



/i + fcI>L = c n n + c s m. 



(18a) 



This equation holds both in the PI and P2 phases. In 
phase P2 no,_ ^ and a second equation emerges from 
the requirement (6_(0)) = 0, which reads as 



/i — fit^L = c n n — c s m. 



(18b) 



Using Eqs. ([IS]). (fT7]) . (fTgaj) and (|18b|) one can get the 
form of the spinor in phase P2, which reads 



2 V c s n c 



(19) 



where n c is the total density of the condensate. 

The obtained phase diagram is depicted in Fig. [T] For 
the critical line one gets 



Bc(T) 

4(o) 



1 



3/2 



(20) 



with T c = (47r 2 n) 2 / 3 fi 2 /(3r(3/2)C(3/2)) 2 / 3 2Mfc B being 
the critical temperature at B = 0. Therefore the shift 
exponent is <p a = 2/3. 



The second step is to express the Green's functions in 
harmony with the equation of state as written above. It 
can be shown that the self-consistent nature of the equa- 
tion of state and the existence of the condensate leads to 
RPA like contributions to the self-energies. To see that 
the procedure leads to a conserving approximation one 
needs to treat correlation functions of the density and 
spin density operators in detail, which goes beyond the 
scope of the present paper and will be published together 
with the dielectric formalism generalized to B ^ [l2j . 

We consider the PI phase first. Since n c _ — in this 
phase the RPA like terms are not present in (we re- 
call that S^j" = anyhow in PI). Consequently Q~^[~ = 
[iu> n — fi _1 £k,-] 1 which after a somewhat lengthy calcu- 
lation leads to ([T5]) with 7 = 2 and 70 = 1 and further- 
more to C(T) ~ T -1 . Further critical exponents have 
the values 77 = rjo = 0, v = 1, Vq = 1/2- They are defined 
as follows: ^^"(0,0) ~ l/k 2 ~ v on the critical line, while 
^^"(0,0) ~ l/fc 2 ~''° at the quantum critical point; the 
correlation length £ - l/|a;-a; c | v and £ - 1/\B-B C (0)\ V ° 
at T > and at T = 0, respectively. The dispersion rela- 
tion is real and parabolic with a gap, that vanishes at the 
phase boundary between PI and P2. The gap exponent 
is vz (T > 0), VqZ (T=0), where the dynamical scaling 
exponent z = 2. Choosing a path towards the quan- 
tum critical point such that [B — B c (0)]/T = const > 0, 
one obtains £ n ~(O,0) - T~ u ° z . The slope of the bor- 
der line of the low temperature region [specified below 
Eq. (fl~5)) ] is proportional to l/c s and is about 200 in the 
units of Fig. [T]by choosing the following realistic param- 
eter values: 



= 1.03 • 10" 50 Jm ; 



and the density (4.6 • 10 20 m 3 ). The correction to the 
leading 1/T singularity has a weaker power law depen- 
dence on 1/T (in the model T" 1 / 2 ). This region, which is 
presumably experimentally accessible terminates at the 
line B = B c (0) along with £ n ~(0,0) - T" 3 / 2 . Note 
that while the power law behaviours are expected to be 
generally valid the exponents along the critical line are 
characteristic for the model. Namely, they are those of 
the spherical model. The reason is that the model can 
be formally derived by introducing a species of spin-1 
bosons and taking the limit a — > 00 while the interaction 
parameters c„, c s are going to zero as \ ja. Furthermore 
B, fi and fx m remain of 0(1). 

In phase P2 the complete self-energy matrix ((7]) has to 
be treated. One obtains for the Green's function 11211 



/ \ rs 

(G_~ = (iu n - fi _1 ek) S rs S ai 

- ftr 1 y/n c , r n c ^ s (C n + rsC s ), (21) 



for r,s € {+, — }. The effective interactions are given 
by C n (k,iu> n ) = c„/[l - 3c n n (0 )(k, iu n )], C s (k,iu> n ) = 
c s /[l — 2c s II(o)(k, iu> n )]. The polarization function II(o) 
is the contribution of the bubble diagram of the free gas, 
since £k, s = in P2 as can be seen from Eqs. (|17[) . (|18ap 
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and (|18b[) and reads as 



II(o) (k, iu n ) 



d 3 q 



"q,s 



(27r) 3 iu n - fi _1 ( e k+q - e q ) ' 



for all values of s. By writing c„,c s instead of C n ,C s 
the Bogoliubov approximation of the Green's functions 
is formally recovered. It is remarkable that the GHP 
theorems (j8|) are fulfilled in this model in such a way 
that simultaneously 



E^(k,0) = E£(k,0)«„* o/J , 



0. 



(23) 



Repeated indices are not summed in this equation. If 
(I23p were true in general then it would mean a gener- 
alization of the results by Nepomnyashchii and Nepom- 
nyashchii [l5[ derived at T = in the case of liquid He- 
lium, i.e. for particles with zero spin. Namely, by ana- 
lyzing the infrared divergences they have shown that the 
anomalous self-energy Ei i _i(O,0) = 0. At nonzero tem- 
perature the expected behaviour is that Ei i _i(k, 0) ~ k 
[111 ] . which is valid here for all the nondiagonal self- 
energies, since II(o)(k, 0) ~ ft -1 . 

The denominator of the Green's function (after ana- 
lytical continuation in frequency) reads as: 



A(k,w) = [lo 2 - H- 2 el] 2 (l ~ 2c s II$)(l - 3c„II$) 
- 2H~ 2 e k [uj 2 - H~ 2 el]n c (c n + c s - ScnCjJffi) 

+ 16fi _4 e k c n c s n Cj+ n Ci _. (24) 



The spectrum of elementary excitations is given by 
A(k, u>)\u=E k -i-t k = 0. Three correlation lengths will be 
used to specify regions in the phase diagram with qual- 
itatively different solutions, namely the thermal wave- 
length A = fi/(2AffceT 1 ) 1 / 2 , the Bogoliubov coherence 
length £b = H/(4:Mc s n c .) 1 / 2 , and the length charac- 
terizing phase fluctuations of the order parameter £' = 
Mk-oT/ (4irH 2 n c _). The three regions of the phase P2 are 
illustrated in Fig [31 For the sake of simplicity we restrict 
the discussion to small k values and assume that c s /c n 
is small, that is the case in the experiment Q ■ Region A 
(A ^> £b,£') : The polarization function can be approxi- 
mated as 



n$(k,«) 



2pe k 



(25) 



where p = r(3/2)C(3/2)/4?r 2 A 3 . The excitation en- 
ergies are linear in wavenumber: oj± — c±k. Here 
c\ = nc n /M — 4c s n + ri- / (nM) and c 2 _ — n c s /M + 
4c s n + n- / (nM), where At T = the 

result reduces to that of Ohmi and Machida obtained us- 
ing the equation of motion method Q. In the present 
approximation a damping appears only for nonzero tem- 
perature, and is exponentially small. In higher order a 
Beliaev type damping can arise (see [ll| for a detailed 
investigation of such a damping in case of a gas of spin- 
less bosons). It is assumed that B is large enough that 



b/bm 



"JC- 
0.2. 



r — n 












FIG. 3: The three regions of the phase P2. Region A and C 
are enlarged for the sake of better visualization. The inset 
shows the order parameter as a function of temperature for 
several values of B. 



the interesting region in P2 keep away from the criti- 
cal line of BEC. Then lu + does not exhibit any remark- 
able change when approaching the P2-P1 boundary and 
therefore only w_ will be followed below. In region B 
(£b ^ £'> A) the polarization function is approximated as 



n$(k )W ) 



ft ? k 



2Mc s £|, oj 



The velocity contains only the condensate densities 

2 4n c ,+n c c s 
c_ = ■ ■ 

n c M 

and the mode has a Landau type damping 



(26) 



(27a) 



7k 



Mc s 
inn 3 



k^Tk. 



(27b) 



Region C (£' ^> £b,\) is the critical region. Here the 
polarization function n^*(k, iui n ) takes the same form 

as in region B. For B, T fixed when k — > the mode be- 
comes overdamped and soft, uj- = — i2fifc/(5M£'). For k 
fixed u_ tends to ek when the critical line is approached. 
The size of these regions depends crucially on c s for the 
other parameters fixed. For realistic values of parameters 
the region A and C are very narrow, the region B is en- 
larged and experimentally accessible. Here the velocity 
depends on the condensate densities; the thermal excita- 
tions represent higher order corrections. More character- 
istically the damping is proportional to the wave number 
and the temperature. These features are known for den- 
sity waves in Bose gases with frozen internal degrees of 
freedom 0, 0j], 0j[ and have been obtained also for spinor 
Bose gas at zero magnetization for the density (and spin 
density) fluctuations [17| • The linear dependence on T is 
even valid for gases in a magnetic trap [H, [H, H3] ■ New 
feature of (|27b[) is the dependence on the strength of the 
magnetic field and magnetization. One expects that the 
leading term remains quadratic in B also for a gas in an 
optical trap. 
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In conclusion, we have investigated the equilibrium and 
dynamic properties of the two Bose-Einstein condensed 
phases of a spin-1 Bosc gas. The richness of properties 
of the two phases from the point of view of a many body 
problem is demonstrated by the GHP theorems. They re- 
flect the basic fact that the number of Goldstone modes 
agrees with the independent phase fluctuations in the 
complex plane exhibited by the order parameter. One of 
them (in P2) is the critical mode for the P2-P1 transi- 
tion, the other will become critical when the temperature 
in PI reaches the line of BEC. The first mode (which co- 



incides with the transverse spin wave only when the mag- 
netization is zero) becomes a quadrupolar spin wave and 
develops a gap in PI. It has been shown that the transi- 
tion between the phases PI and P2 at zero temperature 
belongs to the category of quantum phase transitions in 
which a critical line starts from quantum critical point 
when the temperature is raised. The measurement of the 
shift exponent would be of particular interest. 
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